THE IMPACT OF CONTAMINATED RR LYRAE/GLOBULAR CLUSTER PHOTOMETRY ON 

THE DISTANCE SCALE 
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ABSTRACT 

RR Lyrae variables and the stellar constituents of globular clusters are employed to establish the cosmic 
distance scale and age of the universe. However, photometry for RR Lyrae variables in the globular clusters 
M3, M15, M54, M92, NGC2419, and NGC6441 exhibit a dependence on the clustercentric distance. For 
example, variables and stars positioned near the crowded high-surface brightness cores of the clusters may 
suffer from photometric contamination, which invariably affects a suite of inferred parameters (e.g., distance, 
color excess, absolute magnitude, etc.). The impetus for this study is to mitigate the propagation of system- 
atic uncertainties by increasing awareness of the pernicious impact of contaminated and radial-dependent 
photometry. 



Subject headings: globular clusters: general- 
Lyrae — techniques: photometric 

INTRODUCTION 



-stars: distances, Hertzsprung-Russell and C-M diagrams, variables: RR 



Extragalactic photometry for classical Cepheids in 
M33, M81, M101, and M106 implies that at a fixed 
pulsation period the variables brighten (WVj„) with d e- 



Kennicutt et al 



creas ing galactocentric distance (e.g 
19981 ). A debate persists as to whether the trend is 
tied to a metallicit y gradient, or arises from photomet- 
ric contamination dMaiaess |2010t Maiaess et al.l [20091 
2011b : iKudritzki et al. 20121 and references therein). 
The degeneracy exists because the surface brightness , 
density, and metal abundance ( Andrievskv et al.ll2004 
their Fig. 3) all increase with decreasing galactocentric 
distance. A similar behavior is observed in photomet- 
ric da ta for RR Lyrae varia bles in the globular cluster 
M15 (iMaiaess et al.l 1201 lal ). whereby distance moduli 
established for the variables likewise exhibit a radial de- 
pendence (Fig. [IJ. However, no linear metallicity gra- 
dient exists across the globular cluster, and thus the 
degeneracy plaguing the extragalactic Cepheids is ab- 
sent. Consequently, observations of variables near the 
core of Ml 5 suffer from photometric contam ination (see 
also lYannv et all Il994t IStetson etal] 120051 ) . A depen- 
dence of the photometry on the clustercentric distance 
may complicate conclusions drawn from B ailey (period 



on clustercentric distance (Figs. [1] [2]). It follows that 
the color- magnitude diagrams are compromised (Fig. [3]). 
Clustercentric photometry has a direct impact on the 
establishment of Galactic/extragalactic distances (e.g., 
brightness of the Horizontal Branch), the Hubble con- 
stant (-Ho), the age of th e universe, and cosm ological 
models (e.g., a w ~ 2 x iJn. lMacri fc Riessll2009() . 



ANALYSIS 



2.1. RR LYRAE VARIABLES 



Distance moduli for RR Lyrae variables may be in- 
ferred from a reddenin g-free Wesenheit function (Wyi c )- 
Soszvhski et al. I (l2003h cite the following VI C Wesenheit 
function for LMC RR Lyrae (ab) variables: 



VI, 



(-2.75 ± 0.04) logP + (17.217 ± 0.008) 



The empirical slope of the function is consistent wit h 
that predicted by models (|Di Criscienzo et al. l2004h . 



ampl i tude) diagrams for globular clu sters (ICorwin et al 



20081 ISmith et al.ll201ll and see also lKunder et al.ll2011 



who describe the impact of the Blazhko effect on that 
diagram) . 

In this study, photometry for RR Lyrae variables 
in the globular clusters M3, M54, M92, NGC2419, 
NGC6441, and M15 are shown to exhibit a dependence 



The distance modulus for each RR Lyrae variable was 
evaluated via: 

m — M = 5 log d — 5 
Wvi c - Wvi c ,o = 5 log d - 5 
V - Rvi,{V - h) - alogPo -P = IMi 
V - 2.55(1/ - / c ) + 2.73 log P -/? = Mo (1) 



Rvi c = 2.55 is the ratio of total to selective e xtinc- 
tion (|Di Criscienzo et al.l 12004 iFouaue et al.l 120071 and 
references therein), a is the slope of the Wesenheit 
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Fig. 1. — Relative distance moduli for RR Lyrae vari- 
ables in M15, M3, and NGC2419 exhibit a positional 
dependence. The effect of photometric contamination 
becomes pronounced with decreasing clustercentric dis- 
tance. The relative distance moduli were computed ac- 
cording to Eqn. (fT]). 



function, and /3 is chosen to be arbitrary (relative dis- 
tance) . A mean of two estimates for the Wesenheit slope 
(LMC/IC4499) was adopted (a VIc = -2.73 ± 0.03). 
Only RR Lyrae varia bles sampled beyond th e core of 
IC4499 were analyzed (| Walker fc Nemed[l996h . 

Relative distances were evaluated for RR Lyrae 



2008), M3 



variables disc o vered in M15 (ICorwin et al.1 ._ 
(iBenko et all l2006h. NGC2419 (iDi Criscienzo et al l 
201 l h , M92 (lKopacki|[200l. M54 {Tavden fc Saraiedinil 
20001) . and NGC6441 ( fpritzl et alj|2003l) . The periods 



of RR Lyrae (c) variables pulsating in the hrst over- 
tone were fundamentalized via logPo ~ logP + 0.13, 
where Pq is the fundamental mode period. RA/DEC- 
distance diagrams were constructed for variables in M15, 
M3, and NGC2419 (Fig. [TJ). As expected the distances 
become particularly discrepant as the cluster core is ap- 
proached. The stellar density and surface brightness 
increase markedly as a function of decreasing cluster- 
centric distance. For M15 and NGC2419 a similar trend 
is observed, namely that distances for stars near the 



core appear underestimated. A significant fraction of 
the variables in M15 are contaminated, and numerous 
stars are likewise affected in M3. Four stars in NGC2419 
appear particularly compromised (i.e., A[i < —0.5), and 
less-pronounced contamination persists to near the pe- 
riphery. NGC2419 exhibits a low stellar densitj0 rel- 
ative to the other globular clusters analyzed, however, 
the density of stars per cubic arcminute is large since 
the cluster is distant (d ~ 83 kpc, §2.3j) . 

Fig. [5] features parameters for RR Lyrae variables 
in M3, M92, M54, and NGC6441 plotted as a func- 
tion of clustercentric distance. RR Lyrae variables in 
NGC6441 and M92 display an analogous trend to that 
noted for NGC2419 and M15. Stars near the core ex- 
hibit smaller distances than objects occupying the clus- 
ter periphery. Seven stars within r < 0.5' (NGC6441) 
appear particularly contaminated (Fig. [2]). NGC6441 is 
projected along a line of sight where the field density is 
high (£, b ~ 354, -5°). RR Lyrae variables in M54 dis- 
play a similar trend to those in M3, whereby the moduli 
and reddenings computed appear nearer and larger with 
increasing clustercentric distance, respectively. It is un- 
clear whether effects associated with the edge of CCD 
detectors play a role. Thus photometric contamination 
is not the only source of concern that imposes a cluster- 
centric dependence (further research is necessary) . Stars 
in M3 located r > 3.5' appear closer than those near the 
core (Fig. [2]). A period-reddening relation indicates the 
variables are redder than objects located near the core, 
which are anomalously blue (Fig. [2]). The period-color 
relation does not account for the width of the instability 
strip, thus a symmetric distribution consisting of nega- 
tive reddenings is expected for a population exhibiting 
E(B — V) ~ 0. However, deviations from symmetry are 
observed for stars r < 1' (Fig. [2]). 

Table Q] highlights the difference between the distance 
computed using all RR Lyrae variables from those lo- 
cated near the periphery (/2— fi r ). The offset is largest for 
M15 {Ji-Hr ~ -0.25) and smallest for M92. M3 exhibits 
a positive offset, and the average value for the five clus- 
ters featuring a negative offset is p, — ji r ~ —0.1. Admit- 
tedly, the differences evaluated are subjective and un- 
certain owing partly to the neglect of superposed trends 
(M54, Fig. [2]), and an inability to overcome poor statis- 
tics and assess the extent of the compromised photom- 
etry (NGC6441, M92, Fig. HJ). Nevertheless, the results 
imply a potential ~ 10% systematic uncertainty in age. 
The clustercentric dependence therefore hinders efforts 
to determine reliable parameters for RR Lyrae variables, 
and age estimates for the host clusters. The clusters 
and their RR Lyrae variable constituents are likewise 
employed to establish the distance scale, from which es- 
timates of Ho and the age of the universe stem. Globular 
clusters themselves set a crucial lower limit to the age 



1 Globular Cluster Database, W. E. Harris, McMaster U. 
http: //physwww. monaster . ca/~harris/mwgc . dat 
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Fig. 2. — Relative distance moduli and reddenings for RR Lyrae variables in M3, NGC6441, M92, and M54 exhibit a 
dependence on clustercentric distance. Photometric contamination appears most pronounced near r ~ 0'. The relative 
distance moduli were computed according to Eqn. ([I}. 



of the universe. All potential biases (regardless of their 
relative impact) must be assessed in the pursuit of pre- 
cision cosmology. 

2.2. GLOBULAR CLUSTER CMDS 

Color-magnitude diagrams for globular clusters are 
used to determine their distance (e.g., horizontal branch, 
subdwarfs) and age (isochrone fitting). The radial de- 
pendencies noted for the RR Lyrae variables (Figs. [I] [2]) 
promulgate into the color-magnitude diag rams. Color- 
magn itude diagrams are tabulat ed for M15 (lYanny et all 
19941). M2 [Lee fc Carnev lll999l). M4 (llbata et al.l ll999l 



Mocheiska et alj l2002h . u> Cen (|Rev et all [20041 ). and 
NGC2808 (Fig. [3]). Groups of stars sampling near the 
core and periphery of the clusters are overplotted. The 
severity of the contamination depends on the distance 
and stellar density (stars per cubic arcminute), data- 
processing software, and contribution to the seeing ow- 
ing to the instrumentation (CCD, telescope, wavelength) 
and environment (e.g. . ground or space- based observa- 
tions ), etc. (see also Yannv et al. Il994 : Stetson et al.l 



2005). PSF photometry does not correct for multiple 
stars coincident well within the FWHM of the system. 
The clustercentric effects (Fig. [TJ [2]) impose a degener- 



acy which complicates the analysis of multiple popula- 
tions and potential radial trends which are population- 
dependent (e.g., uj Cen), although for the latter mixing 
over such a large timescale and mass segregation like- 
wise require consideration. The signatures of multiple 
populations may be mimicked by photometric contami- 
nation. 

Photometry acquired for M13 from the Abbey Ridge 
Observatory (ARC@) confirms the trends observed in 
Fig. El For example, as expected photometry near the 
core is artificially brighter and less complete owing to 
the increased sky brightness. The ARO's compara- 
tively small aperture telescope, in tandem with typically 
large local seeing, conspire to create a sizable FWHM 
(~ 2.5 — 3.0"). HST seeing is more than an order of 
magnitude smaller, and for that reason the ARO data 
were acquired to provide an empirical means to evaluate 
contamination effects directly. The trends observed are 
not unique to a given dataset cited, but may be present 
(to an extent) in most globular cluster photometry. 



2 The ARO is a 0.3-m robotic optical facility u sed to detect super- 
novae and conduct variable star research (e.g., [Lane & Gray| [2005l ; 
lLandl200Sft . 
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Fig. 3. — Color-magnitude diagrams for M15, M2, M4, uj Cen, and NGC2808. Black dots define stellar constituents 
closest to the cluster core, whereas red dots characterize a subset of stars beyond the core. The mo rphology displayed 
in the color-magni tude diagrams is sensitive to the clustercentric distance (see also discussion in lYannv et all 11994 



Stetson et al. 2005) 



2.3. 



DISTANCE TO NGC2419 



A distance to NGC2419 may be derived which is 
less sensitive to the trends highlighted in Fig. [TJ Ab- 
solute distance s may be secured by anchoring the 
Soszvhski et al.1 (|2003l ) RR Lyrae VI C Wesenheit func- 
tion to the LMC distance established by lMajaess et al 
(|2011al no, lmc = 18.43 ± 0.030: 



W VIc (RRL,LMC) = (-2.75 ±0.04) log P 
+ (17.217 ±0.008) 



3 A mean computed from ~ 5 X 10 2 LMC distance estimates tab- 
ulated in the NASA/IPAC Extragalactic Database of galaxy Dis- 
tances (NED -D) supports that result jSteer fc M adorc 201l], see 
also Fig. 2 in lFreedman fc Madorefeoiofh 



W V I c ,0 = alogPo +/? 
W VIr - W V i c fl = 5 \ogd - 5 
(-2.75 ± 0.04) log P + (17.217 ± 0.008) - 
(alogPo +p) = no 

ft can be solved by re-arranging the aforementioned 
equation. The absolute Wesenheit function character- 
izing RR Lyrae variables is thus: 

W vu .q ~ (-2.73 ± 0.03) logPo - (1-21 ± 0.03) 



Maiaess et al. ( 2011al ) derived the LMC distance via a 
universal Wesenheit template, which leverages the sta- 
tistical weight of the entire variable star demographic to 
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Fig. 4. — The universal Wesenheit template leverages the statistical weight of the entire variable star demographic 
to constrain distances and pulsation modes for target populations (e.g., NGC2419). SX Phe, RR Lyrae, and Type II 
Cepheid variables in NGC2419 which lie r > 2.5' are shown (gray dots). 



establish precise (< 5%) distances. The template cap- 
italize s upon HST ([Benedict et al.ll2007t ). VLBA, and 
HIP ( van Leeuwen 12007 ) geometric distances for SX 
Phe, S Scuti, RR Lyrae, Type II and classical Cepheid 
variables. Specifically, the templ ate is anchored to 8 
SX Phe and 6 Set variables (HIP, Ivan Leeuwenll2007l) 
4 RR Lyrae variab les (HIP and HST, iBenedict et all 



2002; Ivan Leeuwenl l2007fl. 2 Type II Cepheids (HIP, 



van Leeuwen 2007V and 10 classical Cepheids (HST, 



Benedict et al.l 120071 ). The template may be refined by 



employing new HST parallaxes for RZ Cep, XZ Cyg, 
SU Dra, UV Oct, VY Pyx, and k Pav (jBenedict et al.l 
20111) . However, certain of the aforementioned stars lack 
reliable mean VI C photometry, which is necessary for 
tabulating Wesenheit magnitudes. The desired observa- 
tions are presently being acquired via facilities operated 
by the AAVSO (T. Krajci & A. Henden/BSM). Type 
II Cep heids detected in the less-crowded oute r region of 
M106 (|Macri et al.ll2006l; iMaiaess et al]|2009l ) were like- 
wise incorporated into the Wesenheit templ ate owing to 



the a vailability of a precise maser distance (|Macri et al 
20061 and references therein). The Wesenheit approach 



obviates the need for uncertain reddening corrections, 
and thus the maser distance is used directly. Admittedly, 
the approach supersedes one of the authors (DM) prior 
first-order (yet ultimately incorrect) effort to unify vari- 
ables of th e instability strip in o rder to establish reliable 
distances ( Maiaess et al.ll2011al and references therein). 

Equation [2] may be employed to establish the dis- 
tance to NGC2419, however, a preferred approach is to 
simultaneously use the SX Phe, R R Lyrae, and Type 
II Cepheid variables discovered by iDi Criscienzo et aT 



(|201ll ). An advantage of the universal Wesenheit tem- 
plate is that known pulsation modes for RR Lyrae (ab) 
constituents may be employed to constrain that (poten- 
tially ambiguous) parameter for the SX Phe variables. 
The universal Wesenheit templat e is shown as Fig. 
The template features Galactic (jMaiaess et al] 12011 
and references therein) and M106 data, in addition to 
LMC variables cataloged by OGLE (e.g.. lUdalski et al. 
19991: ISoszvhski et al.l 12008) which were incorporated 



into the temp l ate usi ng the LMC distance derived by 



Maiaess et al.l (|2011al) . Designations for several vari- 
ables tabulated in the OGLE survey are inconsistent 
with their Wesenheit magnitudes (e.g., overtone classi- 
cal Cepheids on the fundamental mode ridge, etc.), yet 
these objects constitute the minority (see also the break 
between 5 Scuti and overtone classical Cepheids). 

To mitigate the effects of contaminated photometry 
stars r > 2.5' (subjective) from the core of NGC2419 
were utilized to deduce the distance. Stars exhibit- 
ing anomalous Wesenheit positions were not included in 
the determination. The offset from the Wesenheit tem- 
plate implies ^ > 19.56 ± 0.05 for NGC2419, which 
is consistent with estima tes tabulated by Harris and 



Di Criscienzo et al.l (|2011l) . The distance is merely eval- 



uated for illustrative purposes since the maximum extent 
of the contamination is unclear (Fig. [T]). The analysis 
reveals that the Type II Cepheid (V18) is either con- 
taminated, a binary system, anomalous, or not a cluster 
member (Fig. [4]). The binary hypothesis is supported by 
the discovery of numerous binary classical and Type II 
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CepheidJjl in the Galaxy and LMC (lEvansI Il99ll Il992 ; 
SzabadoJhooi ISoszynski etaD [2008)). The three SX 



Phe variables lie on the Wesenheit ridge tied to fun- 
damental mode pulsators (Fig. @|. The stars are coin- 
cident with the absolute Wesenheit mag nitude for the 



protot ype of the class (SX Phe, see also iMaiaess et al 
2011ah . SX Phe and other metal-poor population II 6 



Scutis lie toward the short-period extension of the We- 
senheit ridge characterizing population I 5 Scutis. 

2.4. CONCLUSION 

The analysis reaffirms that a (deleterious) cluster- 
centric distance dependence exists in photometry associ- 
ated with the globular clusters M15, M2, M3, M54, M92, 
NGC2419, NGC6441, u Cen, M4, M13, and NGC2808 
(Figs. [U [21 O- Distances inferred from RR Lyrae vari- 
ables and cluster stars are aversely affected. The im- 
pact may be reduced by remaining cog nisant of the 



trends exhibited in Figs . [TJ [3J [3] (see also lYannv et al, 
1994 IStetson et alll2005l ). The principal impetus of the 
present work is to increase awareness concerning clus- 
tercentric photometry in order to mitigate the propa- 
gation of biases into research on RR Lyrae variables, 
and to serve as a catalyst for additional discussion and 
continued research. Photometric contamination, in har- 
mony with the challenges of establishing precise, stan- 
dardized, multi-epoch, multi-band photometry: are fore- 
most among sources of u ncertainty asso c iated with the 
distan ce scale (see also Maiaessl 201fj Maiaess et al 



2011bl ). Admittedly, the radial dependencies in the glob- 
ular cluster photometry were par tly overlooked by one 



2011aL 



of the authors (DM) previously (jMaiaess et al 
and references therein) . 

The conclusions do not mitigate the broader im- 
portance of the studies discussed. For example, the 
datasets contain sizable numbers of newly discovered 
variable stars (e.g., SX Phe, RR Lyrae, and Type II 



Ceph eid variables, IClement et al.l 12001b ISamus et al 
2009)0 The pulsation periods determined place in- 
valuable constraints on models since period-change 
is a proxy for stellar evolution (dR/dt). Such anal- 
yses have hitherto yielded seminal results for RR 



Lyrae and Cepheid variables ( 


Lee & Carnevl 19991: 


Turner et al. 2006t Kunder et al. 


1201 ll Jurcsik et al. 


2012; Ncilson et al. 20121). 
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Table 1 
Relative RR Lyrae Distances 



Cluster 


M - Mr 


M15 


-0.25: 


NGC2419 


-0.11: 


M3 


0.07: 


M54 


-0.08: 


NGC6441 


-0.08: 


M92 


-0.01: 



1 fi — fir defines 
the distance spread 
between the average 
computed using all RR 
Lyrae variables and 
only those near the 
periphery. 
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